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Abstract 

In this paper we show, how a straightforward and natural apphcation of a pair of 
fundamental identities valid for polynomials orthogonal over the unit circle, can be 
used to calculate the determinant of the finite Toeplitz matrix, 



with the Fisher-Hartwig symbol, 

w{z) = C{1 - z f+'f^il - 



w{z) dz 



n—l 



j,k=0 



1, a > -1/2, p e 



Here C is the normalisation constant chosen so that = 2^ • We use the same approach 
to compute a difference equation for expressions related to the determinants of the 
symbol 



w{z) 



a symbol important in the study of random permutations. Finally, we study the 
analogous equations for the symbol 



w{z) 



M 

n 

a=l 



Z- ttn 



9c, 
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1 Introduction 

The large n behavior of the determinant of n x n Toephtz matrices have seen many diverse 
physical applications, from the Ising model jTHj, Random Matrix Theory String Theory 
Pm ITBj and Combinatorics jHl El El- For an arbitrary weight or symbol w{z), the n x n 
Toeplitz matrix is defined by 



For smooth symbols w the asymptotic behavior is a consequence of the classis Szego Limit 
Theorem. For symbols that have jumps or other kinds of singularities, the asymptotic behav- 
ior is harder to determine. The earliest general conjecture for asymptotics of determinants 
for such singular symbols dates back to the formulation of the the Fisher-Hartwig conjecture. 
One of the key ingredients in the proof of the conjecture was the exact computation of the 
determinants for symbols of the form 



These symbols are called symbols with a pure Fisher-Hartwig singularity. The exact com- 
putation was then combined with a localization technique to find more general answers |21. 
A good account of these results can be found in 3J. It is worth noting that in the case 
of only jumps {a = 0) the matrix reduces to a Cauchy matrix and thus the determinants 
are straightforward to calculate. The more general case was done by factoring the Toeplitz 
matrix into a product of triangular and diagonal matrices [71. 

Although, as mentioned above, asymptotic behavior in the case of smooth symbols is 
given by the Strong Szego Limit Theorem, there are applications where more exact in- 
formation is desired. In the analysis of problems involving random permutations certain 
difference equations for quotients of the determinants also arise. These are related to the 
symbol w{z) = e*(^+^/^). For more results about the connection to random permutations and 
random matrices see |Hl El El • 

The purpose of this paper is to re-derive these results and obtain new ones using some 
very simple ideas from the theory of orthogonal polynomials. While the analysis of such 
singular symbols and smooth ones is generally very different the two classes of examples 
share the property that the derivative of their logarithms are rational functions. This is 
what allows us to treat these classes as two examples of a general theory. 

In the course of his investigation into the inversion of finite Toeplitz matrices, Szego 
introduced orthogonal polynomials supported on the unit circle and the corresponding Szego 
kernel. To fix notations for the rest of the paper, let 




(1.1) 





(1.2) 



satisfy the orthogonality condition 




m, n = 0, 1, 2, ... 



(1.3) 
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Note that with the normahsation on the weight, fco = 0o(O) = 1- From the orthogonahty 
condition, we find, jTTj 

knZCpniz) = A;„+l0„+l(2;) - 0n+l(O)0*+i(^) (1.4) 
/^n0n+l(^) = fc„+l#„(2;) + 0n+l(O)0*(2;), (1-5) 

and 

A;„0„(O)0„+i(2;) + A;„_i0„+i(O)0„_i(2;) = (/i;„0„+i(O) + k^+icpniO) z) 0„(z) (1.6) 
by ehminating 0*(-2) from (1.3) and (1.4). The * operation is defined as follows: If 

7r(z] := > ttiZ^ 



j=0 



then 



j=0 



The coefficient of , In can be obtained from the first order difference equation. 



L+l In ^ f 0n(O) \ 0n+l(O) 



kn+l kn \ kn J ^n+l 

if we are able to find fc„ and 0„(O) from the weight. A simple calculation shows that the 
Toeplitz determinant is 

There are, as one may expect, very few cases where explicit formulae were found for the 
Toeplitz determinants. However, it is clear that in order to find information about the 
Toeplitz determinants, a possible approach is to find explicitly the kjS. 

Our first step in that attempt is the derivation of two fundamental identities valid for all 
z E CF'^. We first introduce the functions An{z) and Bn{z) given by 

An{z) = n——+l——Z 0n(O0n 0^ OC?^ 

kn 0n(O) y|^|=i Z-i 

kn An{z) n , f V\z)~v\i) 



Bn{z) = -r^^^ ^ ' ' /"> n(O0n(eV(O^e, 

Kn-\ Z Z J\S,\=\ ^ ~ K 

where w{z) =: exp(— t>(z)). Then we claim 

D / \ I D I \ ^n-\An-\\z) kn 0„_l(O) n — \ , 

Bn{z) + Bn-l{z) = + ———An-l{z) V (z), [Si 

kn-2 Z kn-2 0n(O) Z 
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{Bn+l{z) - Bn{z)){zkn+l(pn{0) + fc„0„+i(O)) 

kn-2 0n(O) 

The equations (Si) and (5*2) are the circular analogues obtained earlier jSJ in the case 
where polynomials are orthogonal with respect to weights supported on the interval [a,b]. 
The functions An{z) and Bn{z) appear in the differentiation formula 

^ + Bn{z)\ = An-l{z)(Pn-l{z) (1.9) 



d , r. /X kn-lAniz) kn 0n-lO \ 

— + Bn-l{z)-- ———An-l{z)](j)n-l{z) 

dz kn-2 Z kn-2 <Pn{0) ^ 

kn 071-1 

{0)An-liz) 



;i.io) 



kn-2 0n(O) Z 

which can be thought of as generalised "creation" and "annihilation" operators. Equations 
(1.9) and (1.10) can be derived jT^j in a manner completely analogous that of the interval 
case. 

The equation (5*1) is simply found by the orthogonality condition and the circular ana- 
logue of the Christoffel-Darboux formula. Let 



(Pniz) 
(Pn-liz) 



Regarding (5*2), this is found by re- writing (1.9) and (1.10), and the recurrence relations as 

= Mn{z)^n{z) 
^n+l{z) = t/„(z)*„(z), 

and demanding that above over-determined systems are compatible; entirely analogues to 
what was done in It is now clear that if w'/w is a rational function then, (5*1) and 
(5*2) will supply the basic equations for the determination of fc„ and 0n(O). Furthermore, by 
eliminating 0„_i(z) from (1.9) and (1.10), we arrive at a second order differential equation 
satisfy by 0„(z) : 



where, 



Y"{z, n) + P{z, n)Y'{z, n) + Q{z, n)Y{z, n) = 0, (1.11) 



Piz,n) = -!^-v'{z)-^ (1.12) 

Z An{z) 

Q{z,n) = BM-Bn{z)^^+Bn{z)Bn-^{z)-'^^^^Bn{z) 

kn 0„-l(O) ( \ n ( \ ^ kn-l(t)n-l{^) An-l{z)An{z) 

-An-l[Z)nn[Z) + T~rn\ • ^--'-'jJ 



kn-2 0n(O) kn-2 0n(O) Z 



We now give one more version of (5*1) and (5*2) which will also prove to be useful in what 
follows. It turns out to be profitable to use the alternative parametrisation, 

0n(O) 



k, 



n+l 
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and the equations (5*1) and (5*2) become, 



Bn+l{z) + Bn{z) = \ V [Z). (Ti) 

m„_i Z TTLn-lSn Z 



S TTi 

{Bn+l{z) - Bn{z)){z + Sn) = m„v4„+i(2;) —An-l{z) - 1. (T2) 



Sn-l mn-2 



2 The Pure Fisher-Hartwig Symbol 

Let us return to the pure Fisher-Hartwig symbol, 

wo^,p{z) = ii-zr^'''{i~i/zr-''' 

where the arguments are chosen so that the first factor is analytic in the interior of the disc 
and one at z = and that the second factor is analytic outside the circle and one at infinity. 
The Toeplitz determinants for this symbol have been previously computed exactly. This 
was done using a hypergeometric approach and a factorization of the corresponding Toeplitz 
operator discovered by Roch. See Theorem 6.20 in [3] for the factorization. An alternative 
direct approach was also done by Widom many years ago, but only recently in printed form 
We now give a difference equation derivation. 

For our purposes we normalize w so that Wq = l/27r. Note also that v' is rational and 
that 

, . , 2a a — ip 

v{z) = - + 

z — I z 

v'{z) - f '(0 _ 2a ^i(3-a 



and thus 



z-i {z-m-i) zi 



An{z) = + {2a - {a + + (2.1) 

kn 0n.(O) Z-l 

B^{z) = ^P-<^-n ^Ml-Ln) ^k^Mzl^ 

Z Z — 1 kn-1 Z 
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where 

In := / Ui)m)'^di 

Now we have 5 unknowns; /„, J„, kn and 0n(O). Notice that An{z) is of the form 
Ci + 02/(2; — 1) and Bn{z) is of the form di/2; + d2/(-2^ — 1)- We now show how five quite simple 
steps yield the difference equation for the knS. 

Step 1. Tcike the limit 2; ^ 00 in (Ti). 

Since Bn{z) has no constant term the constant term of An must be zero. Thus we 
immediately know that 

2amn-l{Jn - In) 



rJz - 1) 



Step 2. Compare the coefficients of I/2; in (Ti). 
A simple computation yields 

2a{Jn-In) 

zp — n — a — 

rn 

and this implies 

mn-i{-a + ip - n) 

This of course also gives information about and at this point we can conclude that 

2n(l - L„) - (\ + i3 - n 



BJz) 



z-1 



Step 3. Compare the coefficients of 1/(2; — 1) in (Ti). 
From this we have 

T 

Aa — 2q;(L„+i + L„) — a-\-il3 — n — 1 = — —i—a -\- ifi — n) -\- 2a. (2.3) 

We will return to this equation in a moment. 

Step 4. Take the limit 2; ^ 00 in {T2). 
This limit shows 

-1 - 2a(L„+i - Ln) = -1 
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or that L„+i = L„. A direct computation shows that Lq — Thus Bn{z) — and 

returning to step three we see that we have an equation for r„, namely, 

—a — id — n — 1 — — —i—a + id — n). 
rn+i 

Step 5. Compare the coefficients of l/{z — 1) in (T2). 

We have that s„ = ^-y^ so this residue produces the equation 

1 H = iTT'ni^ — zp + n + 1) m Aa + zp + n — 1) 

n I n ' n 

and since we have an expression for the quotients of the r„s this becomes 

2q; + 2n + 1 = m^{a + i(5 + n + l){a — i(5 + n + 1) — m^_i{a + i(5 + n){a — i(5 + n), 

(2.4) 



We easily verify that, 



= (n + l)(2a + n + l) 

" (Q; + i/3 + n + l)(Q;-i/3 + n + l)' ^ ' ^ 



solves (2.9) with the initial condistions = and = 1. Prom this it follows 

2 _ r(2a + 1) r(ct + + n + 1) T{a-ip + n + 1) 

" ~ r(a + l)r(a + i/5 + 1) r(n + l) r(n + 2a + l) ^ ' ' 

u _ r(2a + l) |r(^ + ^/j + n)p 

1*^"^ ^' ~ r(n + l)r(n + 2a + l) |r(a + ^/3)|2 ^^"'^ 

= 4=^^.- (2.8) 
n-\- a-\-ip 

Incidentally /„ and J„ can also be easily determined. In Step 1, the vanishing of the constant 
term in A^^z) implies 

2aJn - (a + i/3)/„ 
n H = U, 

and when combined with the first equation in Step 2 gives 

In — "^n 

(a + i/3 - n)rn 



Jn — 



2a 
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3 Toeplitz determinant and discriminant. 

The computation of the Toephtz determinant A„, is now immediate 

n— 1 

j=0 J 

+ l)G{n + 2a + l)G(a + ip + l)G{a ~ ip + 1) 
~ G{2a + l)G{n + a + i(3 + l)G{n + a - ip + 1) 

where G{z) is Barnes G— function. Let 

n 

T^n{z) =-fYl{z - Zj{n)), 

i=i 

where Zj{n) are the n (simple) zeros, then the discriminant is 

\2 



l<j<k<n 



Indeed (see (5.15) of ^) 



\n-l " 



i=i /=i 



n-l 



n— 1 n 



n-l n j^]^ 

The only missing information required for the computation of -D[0n] is 

n 

1=1 

The derivative of 0„(z) is (see (1.10)), 

(f)'^{z) = An{z)(f)n-l{z) - Bn{z)(j)n{z), 

where in our example, 

An{z) = -^-T' a-n = mn-i{i(3 - a - n) 



Now 



1=1 1=1 ' 



0n(l)^ 



from the fact that 4>n{z) = /c„ 11"=! (-^ — zi). To determine 0n(l), evaluate both sides of (3.4) 
at z = 1; keeping in mind that the l.h.s. is regular at every finite z. Therefore the residue of 
the r.h.s. at z = 1 must vanish: 

0.(1) = ^0„-i(l) 



Now 



Therefore, 



j = l J 



D[(t>n] 



-1)V 



0n(l)^ 



(-i)"(-^)/^(./3-a-n)"|^ (n\j) (n^ 



-1) 



n(n+3)/2 



27r 



72 + « — \ " r" ^ 



T(n + 1 + a + i/3) r(a - + 1) V{2a + l)r(r2 + 1) 
T{n + l + a-il3)T{a + il3 + 1) V{n + 2a + l) 



|n + a — 
2^ 



I n— 1 



Note that according to general theory p!Tlj. 

lim A;„ = K > 0, 

n— »oo 

and in particular 

lim 0„(O) = 0. 

n— >oo 

Indeed this is the case for the Fisher-Hartwig symbol. 



(3.5) 



(3.6) 



(3.7) 
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To determine as n — > oo, we give here some preliminary results: 



lim k„ — 



VT{2a + 1) 
\r{a + i/3 + l)\ 



In— 1 



r(a - i(3) 



n 



A/r(2Q; + l) 



G{2a+l) 



Therefore, 



\V{a + il3 + l)\ G{2a + l) 



r(2a + 1) 



\a + i(5\ + + \C\a + i(5\\T{a + i(3)\'^ 

4 Differential Equation. 

From the general theory (pn^z) satisfies 

Y"{z) + P{z, n)Y\z) + Q{z, n)Y{z) = 0, 

where 

„, , 1 — n — a + iB 2a + 1 
P{z,n) = - + 



n 



(3.8) 



Q{z,n) 



z z — 1 ^ 

n{a + i(3 + l) 



z(z - 1) 

The general solution of the differential equation is 

A 2Fi{—n^ a + i[3 + 1-^1 — n — a + iP; z) 
+ B z^'^''-''^ 2Fi{n + 2a+l,a-ip;n + a-ip + l;z), 

and since we know that 0„ is a polynomial we have that 

<Pn{z) = A 2Fi{—n, a + ip + 1; 1 — n — a + iP; z) 

where 



(4.1) 



(4.2) 



T{2a + l)r(a + i/3 + n + l)T{a - i(3 + n + 1) T{n + a- if3)T{a + il3 + 1) 
r{a -i(3 + a)r{a + ip + l)r(n + l)r(n + 2a + 1) r{n + a + ip + l)r(a - i(3) ' 
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This value comes from computing the coefficient for the nth term of the Hypergeometric 
function and also using the value of kn given in (2.6). Indeed this particular Hypergeometric 
function was found by Askey to be orthogonal with respect to the pure Fisher-Hartwig 
symbol in a commentary on Szego's collected papers [T]. 

5 The weight Cexp[t{z + l/z)/2]. 



For this example 



t t 

2 ^ 2^2 
v'{z) - v'iCl _ 1 t It 

z-^ ~ ~^2^ ~ ?2^ 



We find, 



z 

jj , . Ln + bn t/2 

Bn{z) = \ 5 

z z^ 



2r„ 



2r„ Ji^i^i <2 



l€l=l 



2 



There are 5 unknowns: r„, m-n, On, and L„,. Again we use the same basic steps as was 
done in the previous example. 

Step 1. Take the limit as 2; ^ oo in (Ti). 

bn = -n-^Sn. (5.1) 
Step 2. Compare the coefficients of 1/z in (Ti). 

t ( '^CLnlt\ 

L„+i + Ln - -\ Sn+i H = n + 1. (5.2) 



2 . 

Step 3. Compare the coefficients of l/z"^ in (Ti). 

t 

an = 2- 



(5.3) 
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Step 4. Tcike the limit as z — > oo in {T2). This hmit yields after simplifying 
L„+i - - ml) - - - m^_J^ = 

which implies, 

where a is an "integration" constant. To determine a put n — and note that m_i = and 
•50 — t^i/t^o — Ti — —Ii{t)/ Io{t), where Ij{t) is the /— Bessel function of order j. So, 

a = Lq- -ri. 



But 



Jo 

nthit) th{t) 



= -7rCth{t) 



2nIo{t) 2Io{t)' 



So q; = 0, and 



Ln^^-Snil-ml_,). (5.4) 



Step 5. Compare the coefficients of in (T2). Again after simphfying 



which implies. 



+ + ) =1, 



Ln - -{sn-\ — ) ^ (3 + n. 

2 V Sn-1 



It turns out that the "integration" constant /3 is also 0. Therefore, 



L„^n+-{sn + ^^) . (5.5) 



2 \ Sn-l 

Step 6. Compare coefficients of in (T2) This gives = 
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Eliminating L„ from (5.4) and (5.5) gives, 

- ^ = <-i (sn + . (5.6) 

Now substitute (5.5), (5.4) into (5.2) and note that a„ = t/2, we get after some simplification. 

Now since s„ = r„_|_i/r„, we see that 

7 = constant. 



2 9 
' n+1 ' n 



I -ml 1 - m^_i 
which implies, 

K-i = ^-rl (5.8) 

where an easy computation shows that 7 = 1. This last equation hold in general jU]. 
Using this on (5.6) gives the discrete Painleve II pTl ITB] 

'in Tn 

r„+i + r„_i = — — -. (5.9) 

t l-rl 

6 Another non-linear difference equation for r^. 

In [inj the authors considered a symbol of the form 
«;(z) =e*^n ( ' Z1^" = ^>0' -l<-2;«<0, J]^?a^a = -1. (6.1) 

Q=l ^ ' 0=1 a=l 

We have, 

v{z) = -tz-^ ga (ln(z - Za) - In z) 

a 



z-^ z^ ^ z - z^ Z^ 



An{z) = (n + a„)m„_i + m„_i ^ ^ 



bJa) 



Z - Zn 



13 



r,. 



Aniz) n + g Y^5'a + c„(tt) 



m , y r ^ — ^ 



rrin-iz z — ' — -^o 



J\^\=l K-Za 



Step 1. ^ ^ oo in (Ti) 



a„ + n = -^^^^i = (6.2) 



Step 2. Residues at z = of (Ti) : 



On = ^ + 2^ ■ (6-3) 



Step 3. Residues z — Za of (Ti) : 

S-a + C„+i(Q!) + C„(Q!) H = . (6.4) 

Za 

Step 4. 2; — i> cxD in (T2) : 
1 + a„+i - a„ + (cn+i(Q;) - £„(«)) = ml{n + 1 + a„+i) —ml_i(n - 1 + a„_i) 

which is simphfies to: 

trlsn - trl^i + (cn+i(a) - c„(a)) = 

using (6.2) and m'^_i = 1 — r^. "Integrating" the above gives: 

trnr„+i - ^ c„(a) = A(a), (6.5) 

a 

where X{a) is the n— independent integration constant and depends on t and the set of z'^s. 
Step 5. Residues at 2; = of (T2) : 

Ebn+i{a) - bn{a) 
) 

a 
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and upon "integration" gives, 



a 



the same as (6.3). 

Step 6. Residues aX z = Za of (T2) : 

{Za + Sn) (cn+l{a) - Cn{a) + ^""^^^'^'^ J ^ 77^2 f!l_^2^_^^^__^ (ct) . (6.7) 

Let's see what all this means when there is just one Za- 
If we use Step 3, along with (6.5) and (6.6) we have that 

g + tr„+2r„+i + A + tr„r„+i + \- g- Sn+it -n-1 = —{-g - sj - n), 

or, after simplifying, 

trn+lTn + +trnrn+l + 2A - Sn+lt - U - 1 ^ 1/ Sn - tZa - nZa/Sn, 

A httle more algebra yields 

_ {nrn - t)za - 2A + (n + 1) 

fn+2 -rfn— -7- 2 \ 

an equation very much like the one found for the previous example. 
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